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ABSTRACT 

In this paper, we present a simple, yet versatile, analytical model of one-dimensional photonic crystal (ID PC). 
In our theoretical model, we take into account direction of propagation and therefore do not neglect anisotropic 
nature of photonic crystals. We derive analytical expressions for mode spectrum and density of states in ID 
I photonic crystal. With those formulas, we obtain mode spectrum characteristics, which depict formation of 
7-H photonic band gap and reveal properties of photonic crystals. 
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^ ■ 1. INTRODUCTION 

' Photonic crystals (PCs) , proposed by Yablonovitclf^ and John,'2l represent a novel class of optical materials which 
^> . allow to control the flow of electromagnetic radiation or to modify light-matter interaction. These artificial 
structures are characterized by one-, two- or three-dimensional arrangements of dielectric material which lead 
to the formation of an energy band structure for electromagnetic waves propagating in them. One of the most 
interesting features of photonic crystals is associated with the fact that PCs may exhibit frequency ranges over 
which ordinary linear propagation is forbidden, irrespective of direction. These photonic band gaps (PBGs) lend 
themselves to numerous diversified applications (in linear, nonlinear and quantum optics). For instance, PBG 
structures with line defects can be used for guiding light. Similarly, as it has been predicted and confirmed 
f*) i' experimentally, photonic crystals allow to modify spontaneous emission rate due to the modification of density 
O . of quantum states. In particular, it is well known that the density of states grows at the edge of the photonic 
c/3 ' band gap of the PC. This allows us to predict a higher optical gain, but on the other hand a higher level of noise 
in light generated in PC-lasers operated at a frequency near the band gap. 

It is our purpose here to present a simple, yet versatile, analytical model of one-dimensional photonic crystal. 
■ The investigated photonic crystal consists of N elementary cells with alternating layers of different indices of 
CIh. refraction and different widths. Presented model takes into consideration direction of propagation and therefore 
does not neglect anisotropic nature of photonic crystals. We derive analytical expressions for mode spectrum 
, and density of states that allow to analyze properties of one-dimensional photonic crystals, especially inhibition 
^ ' or amplification of spontaneous emission due to modification of density of states. The described model allows 
'n]" to examine basic physics behind PBG structures and as an analytical model - can provide an excellent starting 
point in numerical calculations for much more complicated structures. 
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^ ; 2. EFFECTIVE RESONANCE APPROACH TO PHOTONIC CRYSTALS 

' The simplest case of a photonic crystal is one-dimensional photonic crystal. It is a structure built of alternating 
' layers of linear, uniform and isotropic materials of two different refractive indices ni and 712. Its section in xz 
. . plane is shown on figure [1] The structure has rotational symmetry around z axis, all layers extend towards 
^ infinity in x and y directions (direction y is perpendicular to section shown on figure [T]). We will consider 
a photonic crystal built of finite number of layers. A primitive cell in this case is a pair of subsequent layers 
with refractive indices ni and rt2. The number of primitive cells "right" from z — plane will be and to the 
Ci I "left" will be , where directions "right" and "left" are in accordance with those on figure [1] We will introduce 
primitive cell index m as depicted on the figure: layers with z G [0, A] constitute primitive cell with to = 0, and 
the next (along with the z axis) pair has m = 1, the first one has to = — 1, the first primitive cell "left" from 
z = plane has to = — 1 and the index of the last one is rn = —N^. We also assume (for simplicity), that the 
crystal is surrounded by a material with refractive index ni. 

The idea of the described method'^ is presented on figure [51 It aims to determine effective reflection from 
a part of photonic crystal on a layer's boundary and "replace" it with a mirror of the same reflection. Thus, 



00 



C/2 



o 

■ 1— I 



X 



Figure 1. Structure of one-dimensional photonic crystal (section). 



we can say, that chosen crystal's layer in which the density of states is investigated, is inside an effective FP 
cavity. This method can describe properties of finite photonic crystals, and allows to examine changes of these 
properties in different layers of the crystal. In general case of photonic crystal, precise description of density of 
states without approximations requires difficult and time consuming calculations. However, for the simplest case 
of one-dimensional photonic crystal it is possible to follow the derivation given below in the further part of the 
paper. 

In our model we assume that in each of layers of the photonic crystal Maxwell equation^ have solutions in 
form of plane- waves: 

E = Ene'''^-"^' (1) 
with wave vector k bound with angular frequency w by dispersion relation 



k' = 



(2) 



where c is speed of light in vacuum and n is refractive index of the material of the layer. We can relate fields in 
different layers using continuity conditions, what is especially simple in case of one-dimensional photonic crystal.l^ 
In every ni layer we represent electric field as a superposition of coupled plane-waves: 



and similarly in every n2 layer: 

Em,2^e = Cra,,e'''-''+"^''-"^' + d„,,e*'=="^-*'3^^~'"'* , z G [ mA + a, (to + 1) A ] 



(3) 



(4) 



where e denotes polarization TE (electric field perpendicular to plane of incidence) or TAl (electric field in plane 
of incidence) and 



kl , J = 1,2 



(5) 



Continuity conditions lead to equation relating amplitudes in consecutive primitive cells: 
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Figure 2. A layer of one- dimensional photonic crystal as a resonator. 

where Mm,e is translation matrix of mth primitive cell. It is an obvious conclusion, that if amplitudes of plane- 
waves outside the photonic crystal are indexed with m = Nj^ for "right" and m = —N^ — 1 for "left", then 



, + 1,, ...M_2,.M_i,,)" 



(7) 
(8) 



with the amplitudes outside the crystal bound by external conditions. ([7]) and are general relations, that 
allow to perform calculations for defected structures. For a perfect photonic crystal we have 
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and 
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As a result, ([7]) is simplified to 
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(for more elegance the index R is omitted). Matrix is unimodular, therefore is obeys the Chebyshev identity: 

(20) 
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r/e = arccos- 



(21) 



what allows to simplify even more. Relation (or ([7]) and ([5]) in more general case) is a set of two linear 
equations.SI Amplitudes aN,e and bN^e are bound by external conditions, in this case it is just 







(22) 



because we assume excitation in the photonic crystal (or at least on the other side of it) and there is no reflection 
in the infinity. Because ao,e is the amplitude of the wave incident on layer boundary in z = a, 60, e can be treated 
as refiected wave. Therefore, we can calculate reflection coefficient of the whole part of photonic crystal reflecting 
ao,e wave, it will be 

bo ee-*^i" 



Knowing p^ . ([^0]) and (1^^ (and again marking index R for "right" part of the crystal), we find that 
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Because the structure is symmetrical, we instantly have expression for "left" part reflection coefficient 
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(in general case of defected crystal it is not possible to obtain such a concise expression for at least one reflection 
coefficient). 

Having calculated reflection coefficients, we treat the layer of photonic crystal as a resonator, with mirrors 
constituted by the remaining layers of the crystal and described by derived formulas. 



3. MODE SPECTRUM AND DENSITY OF STATES IN A LAYER OF 
ONE-DIMENSIONAL PHOTONIC CRYSTAL 

Let operator Q relate electromagnetic plane-wave E'^^ exciting the resonator to effective plane-wave E'^f^ with 
the same wave vector, that appears in the resonator due to excitement: 

E^ff = QE^"" (26) 

We call a plane- wave E^ [k^ with polarization e and wave vector k a mode of the resonator if it satisfies eigenequa- 
tion of Q: 

QE,(Pj ^ K,(Pj E,(Pj (27) 

If the resonator is excited with plane- wave E^ (^tj , effective distribution of electromagnetic field is superposition 

of plane-wave A^^fc^ E^(^}ij with its coupled mode (i.e. plane wave with wave vector [ /c^^, fc^, — fcz ]), which 
appears because of reflections. Therefore, introduced definition is equivalent to "classical" definition of a mode. 

We introduce mode spectrum ^fc^ , related to distribution of modes in wave vector domain. Thus, it can 
be used to calculate the number of modes N{k) with propagation constants less than k: 

N{k) = V / dk' k'^ [ dfl' pe (k') (28) 

Differential of N{k) is known as density of states. From (pS)) it can be found, that 



dN{k)dk ^ k^dkJ2 dnp,(tj 



(29) 



Because operator Q can be used to find distribution of modes, mode spectrum (je^ has to be proportional to 
its eigenvalues: 

pJk)^NpAjk) (30) 



From (P^l we can calculate, that in free space, where Ag{^kj ~ 1, 

dNfs{k) dk = k^dk ^ / dnNp:^ SjrNpk^dk (31) 

g Jatt 

It is a well known fact, that 

1.2 

dNfs{k)dk^ —dk (32) 
Comparing (I3ip and p2p we find the normalization constant Np: 

Np = (33) 

Finally, mode spectrum is defined as 



{2^Y 

Mode spectrum depicts influence of a FP cavity on a plane-wave. It also contains information on anisotropy of 
the structure, and therefore is a magnitude more fundamental than density of states. 

To find mode spectrum of a ni layer of one-dimensional photonic crystal we need to find operator Q for a FP 
resonator filled with medium with refractive index rii, built of two parallel mirrors at distance a from each other. 



with reflection coefficients g and r^^g. If the resonator is excited with a plane- wave with wave vector k and 
amplitude of electric field Eq (this is a scalar problem), then amplitude of electric field of effective plane- wave 
with wave vector k can be found by summation of two infinite series 

oo oo 

E^ff ^Eo + Y, rR,,rL,ee^'^''Eo + ^ 4^,r2_,e-2^^"i?o (35) 

where /3 = /3i defined by ([S]) and * means complex conjugation. Comparing to definition (pS)) we can conclude, 
that the operator Q for the resonator is: 

Q = ^-\rR.r,,\' (36) 

1 + \rR,erLA - 2Re{rij,,rL,,e2'/3«} 

and from ([M)) we obtain the formula for mode spectrum of the photonic crystal's layer: 

f fc) = A (37) 

To find density of states we can use definition (P^)) . In spherical coordinates (3 — kcos-d, therefore 
dNik) = e Cm' Bind' ^dip' y A 1 - \rR,erL,e\^ 

Integration over ip' can be done, because of rotational symmetry of the crystal. Additionally substituting 
u = cosi?' we reach the following formula: 

dN{k) ^Jf- f duY ^-\rR.^rLA^ 

The above expression is true for any resonator with two parallel mirrors. In particular, taking tr^^ = r^^^ = 
and rii = 1 we obtain expression for density of states in free space: 

dNfs{k) = ^dk (40) 
If both r^^g and vr^^ are constant, mode spectrum takes form of regularly distributed "peaks" (see figure [21). 



4. VARIATION OF MODE SPECTRUM IN A ONE-DIMENSIONAL PHOTONIC 

CRYSTAL 

Presented method of mode spectrum calculation allows to analyze its space dependence in photonic crystal and 
therefore also the crystals properties. For example, it is possible to examine the photonic band gap edge in 
different layers of the crystal (figure U]). It can be clearly seen, that in layers (sufficiently) close to the end of the 
photonic crystal full photonic band gap does not form. The reason of this behavior is that reflection coefficients 
and rp depend on number of layers surrounding selected layer. If this number is high, amplitude of reflection 
coefficient can be close to 1, while for low number of layers amplitude of reflection coefficient is signiflcantly 
lower. In the second case, quality of resonator is low and its mode spectrum tends to be more similar to mode 
spectrum of free space. Therefore, formation of distinctive photonic band gap does not occur. This regularity 
can be as well observed on flgure [SJ which depicts mode spectrum (for perpendicularly incident plane- wave) in 
different layers of the crystal. Amplitude of mode spectrum at the edge of photonic band gap reaches highest 
values in the middle of the crystal, where amplitudes of reflection coefficients and rp, along with effective 
resonator quality, have the highest values. Another interesting observation, is that photonic band gap in general 
is not symmetrical - amplitudes of mode spectrum at the edges of the gap are not equal. 
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Figure 3. Modes of FP resonator with constant reflection coefficients of mirrors. 
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Figure 4. Edge of photonic band gap in different layers of a ID photonic crystal. 
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Figure 5. Mode spectrum in a ID photonic crystal. 
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